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Symmetry of the remanent state flux distribution in superconducting thin strips: 

Probing the critical state 
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The critical-state in a thin strip of YBa2Cu307_5 is studied by magneto-optical imaging. The 
distribution of magnetic flux density is shown to have a specific symmetry in the remanent state 
after a large applied field. The symmetry was predicted [Phys. Rev. Lett. 82, 2947 (1999)] for 
any jc(B), and is therefore suggested as a simple tool to verify the applicability of the critical-state 
model. At large temperatures we find deviations from this symmetry, which demonstrates departure 
from the critical-state behavior. The observed deviations can be attributed to an explicit coordinate 
dependence of j c since both a surface barrier, and flux creep would break the symmetry in a different 
way. 



I. INTRODUCTION 

During the last years much attention has been paid to 
studies of the magnetic behavior of thin superconductors 
placed in a perpendicular magnetic field, the so called 
perpendicular geometry. On one hand, numerous papers 
have investigated theoretically the critical state of thin 
superconductors of regular shapes like long strips and cir- 
cular disks.ErQ On the other, magnetic characterization 
by measuring the spatial distribution of flux density at 
the surface of superconductors has become quite com- 
mon and powerful. This progress has been facilitated 
by the development of spatially-resolved techniques, such 
as magneto-optical imaging, Hall- probe arrays etc., see 
Ref. ^ for a review. 

In spite of these massive efforts, the task of making 
proper interpretations of a measured flux density distri- 
bution, -B(r), is still one with major difficulties. In par- 
ticular, we are not aware of any simple decisive method 
to judge whether an observed B(r) is consistent with the 
critical-state model (CSM) or not. One could here expect 
that fitting an observed B(r) by profiles predicted from 
the CSM with some .B-dependent critical current density, 
j c (B), would be a straightforward procedure. However, 
this is not so since in the perpendicular geometry explicit 
expressions for the flux density distribution are available 
only for the Bearumodel, j c = const, for a thin striprm 
and a thin disJptl, as well as for a strip with a special 
kind of j c (B)B For a thin strip or disk with a general 
-B-dependent j c , the flux distribution can be calculate' 
only numerically by solving a set of integral equations!]' 

In a recent workj!3 we considered the critical-state 
magnetic behavior of a thin superconducting strip with 
a general independence of j c . It was shown that the 
central peak in large-field magnetization hysteresis loops 
of such samples always occurs at the remanent state, 
B a = 0. An intermediate result of that derivation is 
the prediction of a special symmetry of the remanent- 
state flux density distribution. Since the symmetry is 



independent of the particular j c (B), it may serve as a 
conclusive and easily implementable test for the applica- 
bility of the critical-state model in a given experiment. In 
the present paper we demonstrate using magneto-optical 
imaging how this symmetry in the flux density profile 
can be revealed, and used to verify the applicability of 
the CSM. 



II. SYMMETRY OF FLUX DISTRIBUTION 

Consider a long thin superconducting strip with edges 
located at x = ±w, the y-axis pointing along the strip, 
and the z-axis normal to the strip plane, see Fig. d The 
magnetic field, B a , is applied along the z-axis, so screen- 
ing currents are flowing in the y-direction. Throughout 
the paper B denotes the z-componcnt of magnetic induc- 
tion in the strip plane. The sheet current is defined as 
J(x) = J j{x, z) dz, where j(x,z) is the current density 
and the integration is performed over the strip thickness, 
d< w. 




FIG. 1. Superconducting strip in a perpendicular applied 
magnetic field. 

From the Biot-SavartJ.aw for the strip geometry, the 
flux density is given byErB 



B(x) - B a 



1^0 

2ir 



J(u) du 



(1) 
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Assume that the strip is in the remanent state after a 
very large field was applied. Then, everywhere in the 
strip the current density is equal to the critical value, 
J(x) — sgn(x) J c [B(x)], and the flux density distribution 
satisfies the following integral equation, 



B(x) = 



J c [B(u)} 



; du . 



(2) 



By changing here the integration variable from u to 
\Jw 2 — u 2 , and also replacing x by \/w 2 — x 2 , one ob- 
tains 



B(V^) = >«r J <Wf r r ?) K d u. (3) 



7T Jq X 2 — U 2 



Clearly, this equation for —B(y w 2 — x 2 ) is equivalent to 
Eq. (||) for B(x) if J c depends only on the absolute value 
of the magnetic induction, J c (|i?|). Thus, we conclude 
that 



B(x) = -B(Vw 2 -x 2 ) . 



(4) 



This symmetry of the flux density distribution in the 
fully-penetrated remanent state is valid for any J C (\B\) 
dependence. In particular, it also holds for the Bean 
model, where one finds by simple integration of Eq. (|^) 
that 



B(x) 



2ir n 



(5) 



It follows from the general Eq. (^) that the flux density 
is always zero at x — w/V2. At this point B{x) changes 
sign from positive in the central part of the strip to neg- 
ative near the edges. 

This special symmetry of the flux density profile across 
a thin strip has a trivial analog for the case of a long 
sample in a parallel field. There, the gradient in B is 
a function of the local value of the field. Hence, the 
flux distribution is always symmetric around the point 
xq where the flux density is zero, B(xq) = 0. Hence, one 
can write B{x) = —B(2xq—x) as long as x and 2xq—x are 
within the superconductor. In the perpendicular geom- 
etry the relationship between B and j is non-local, and 
the B-profiles deviate strongly from those in the interior 
of long slabs and cylinders. 



III. EXPERIMENT 

A 200 nm thick film of YBa2Cu307_,5 was gro\pa epi- 
taxially on an MgO substrate using laser ablation.Eil The 
sample was patterned by chemical etching into a strip 
shape of width 2w — 2.5 mm. For the measurements of 
flux density profiles we chose a region free of any defects 
visible by our magneto-optical (MO) imaging system. 




FIG. 2. Flux density profile across the YBa2Cu307_5 strip 
in the remanent state at 42 K. The profile is derived from the 
magneto-optical image shown above. The bright areas on the 
image correspond to large |B|. 



42 K 




FIG. 3. Profiles of the flux density across one half of the 
YBa2Cu307_<5 strip at different temperatures. The strip is in 
the remanent state after a large field had been applied. 
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The imaging system is based on the Faraday rotation 
of polarized light illuminating an MO-active indicator 
film that we mount directly on top of the superconduc- 
tor's surface. The rotated Faraday angle varies locally 
with the value of the perpendicular magnetic field, and 
through a crossed analyzer in an optical microscope one 
can directly visualize and quantify the held distribution 
across the covered sample area. As Faraday-active indi- 
cator we use a Bi-doned yttrium iron garnet film with in- 
plane magnetization.llj The indicator film was deposited 
to a thickness of 5 /zm by liquid phase epitaxy on a 
gadolinium gallium garnet substrate. A thin layer of alu- 
minum is evaporated onto the film allowing incident light 
to be reflected, thus providing double Faraday rotation of 
the light beam. The images were recorded with an 8-bit 
Kodak DCS 420 CCD camera and transferred to a com- 
puter for processing. The conversion of gray levels into 
magnetic field values is based on a careful calibration, see 
Rcf. |l3|. The MO imaging at low temperatures was per- 
formed by mounting the superconductor/MO-indicator 
on the cold finger of a continuous He-flow cryostat with 
an optical window (Microstat, Oxford). 

MO images were taken in the remanent state after ap- 
plying a large field at temperatures of 42, 58 and 82 K. 
An MO image at 42 K and the corresponding flux den- 
sity profile are shown in Fig. ^. The fact that maximum 
trapped flux density is observed in the center of the strip, 
implies that the applied field had been raised to a suffi- 
ciently large value. Furthermore, one sees from the fig- 
ure that the return field of the trapped flux penetrates 
regions near the strip edges, x = ±w, where the field 
has negative peaks. While the intensity in the MO im- 
age does not immediately discriminate between the two 
field polarities, it is readily accounted for by locating 
the boundary where B = 0, also called the annihilation 
zone. One sees also that the flux distribution in the left 
and right halves of the strip are mirror images of each 
other, and therefore we focus only on one half of the 
strip, < x < w. The flux density distribution at higher 
temperatures are shown in Fig. ^. Due to reduced flux 
pinning with increasing temperatures the magnitude of 
the trapped field is reduced substantially. In addition, we 
also see changes in the shape of the flux profile. The spa- 
tial resolution of the method is limited by the thickness 
of the MO indicator film. Therefore, a few data points 
in 5 /im vicinity of the peaks at the strip center and at 
the edge have been removed from the following analysis. 

To test the symmetry property, expressed in Eq. ([!]) , of 
the measured flux profiles we plot in Fig. ^ the absolute 
of the flux density as a function of the new coordinate x' , 

/ _ f x, 0<x<w/V~2 

\ \fw 2 — x 2 , w/V2 < x < w . 

If Eq. ([I|) holds, there should be full overlap of the two 
branches where the measured B is positive and negative. 
One sees from the figure that the overlap is almost com- 
plete for the data taken at 42 and 58 K except for small 



deviations at large \B\. At 82 K, however, there is a 
significant splitting of the two branches over the whole 
range. One may therefore conclude that there is a sys- 
tematic deviation from the CSM behavior at this high 
temperature. 
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FIG. 4. Profiles of the flux density from Fig. |3| replotted 
with new coordinate x' defined by Eq. (0). Open and solid 
symbols correspond to x > w/y/2, and x < w/y/2, respec- 
tively. A large splitting of the two branches at 82 K indicates 
a deviation from the CSM. 

The results from the symmetry analysis of £?-profiles 
are now compared with direct evaluation of j in regions 
with equal \B\. For this purpose the current density dis- 
tributions were calculated from the measured B-profiles 
by the inversion-scheme proposed in Ref. |l3| and devel- 
oped elsewhere." The latter procedure is, in general, 
much more complicated than the analysis of B-profiles 
and it requires knowledge of ^-distributions across rather 
large regions also outside the strip. Furthermore, the in- 
version procedure involves filtering of short-wavelength 
noise in experimental data. 

The current profile at 42 K is shown as a solid line in 
the inset of Fig. |^. One can see an enhancement in j 
in the region near B = 0. In the main figure the data 
are replotted as j versus \B\. Again we see that the two 
branches corresponding to positive and negative B col- 
lapse, which proves existence of the critical state in the 
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strip. The unified curve characterizes the j c (B) depen- 
dence. Also in this plot, like in Fig. [|, there are small 
deviations from data collapse at the largest fields. We 
thus see that the analysis based on j-inversion gives sim- 
ilar results as the symmetry analysis of B profiles. 
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FIG. 5. Sheet current versus the absolute value of the local 
flux density at 42 K. The data are obtained by replotting the 
B(x) and j(x) profiles shown in the inset. A collapse of the 
two branches, which correspond to positive and negative B 
(solid and open symbols), demonstrates that the critical-state 
is established in the strip. The dashed line shows a fitted Kim 
model j c (B) oc (1 + B/Bq)' 1 with B = 150 mT. 



IV. DISCUSSION 



The symmetry condition, Eq. (Q), was derived within 
the framework of the CSM, and is not expected to hold 
when one or more of the basic assumptions of this very 
successful model is violated. The most probable reasons 
for lack of this symmetry are therefore the following; 

(i) Presence of a surface barrier for vortex entry and 
exit. There are vast experimental observations of such 
a barrier in high-T c superconductor crystals. The bar- 
rier leads to an excessive current density in the vicinity 
of edges, which then destroys the symmetry In particu- 
lar, it shifts the point xq where B = towards the strip 
edge. For the geometrical barrier which arises from a 
rectangular shape of the cross sectional area, JJie exces- 
sive current density is of the order of H c i/d£3 For the 
Bean-Livingston barrier it should be larger and scalejsath 
temperature as the thermodynamic critical field, H c £3 In 
both cases, the temperature dependence can differ from 
that of the bulk pinning j c suggesting that deviation from 
the symmetry can be temperature dependent. Moreover, 
it is known from experimentO that the surface barrier 



dominates at the higher temperatures, while bulk pin- 
ning is more important at low temperatures. 

(ii) Thermally- activated creep of vortices leading to a 
slow time relaxation of the flux distribution. The flux 
creep problem for a thin strip with a i?-independent 
E(j) law has been considered in Refs. Under 
constant applied magnetic field the space and time de- 
pendence of the electric field is shown to decouple as 
E(x, t) — f(x)g(t), where f(x) can be found numerically. 
It is also argued that during relaxation of E(x, t), starting 
from some initial E(x,0), the electric field will approach 
the profile given by f(x). From B = —dE/dx it follows 
that if an initial remanent flux profile B(x) crossed zero 
at xq — then during relaxation xo will shift to- 
wards the point where f(x) has the maximum. For the 
voltage-current law E = E c (j / j c ) n , the maximum in / is 
always located at Xq > w/\/2, namely at xp = 0.735w for 
n = 1 and approaches w/V2 as n — > oo S3 This means 
that at smaller n, i. e., at larger temperatures the devia=. 
tions from the symmetry due to relaxation are stronger £3 

Thus, both a surface barrier as well as flux creep pre- 
dict a stronger deviation from the symmetry in B{x) at 
higher temperatures. However, neither of them can ex- 
plain the deviation found in our experiment. Indeed, 
while the CSM predicts that in the remanent state xq = 
w I V2, both surface barrier and flux creep lead to larger 
xo. Such a shift of xq would result in the negative- B 
branch being below the positive- B branch in the |_B|(x') 
plot. However, that is just the opposite to what is shown 
in Fig. f| 

(iii) The last possible reason for the deviation from 
the symmetry is inhomogeneity of the strip, which leads 
to an explicit coordinate dependence of the critical cur- 
rent density, j c {x, B{x)). It may be caused, e. g., by a 
nonuniform chemical composition£j The kind of devia- 
tion shown in Fig. ^ can be explained by a suppressed 
j c near the strip edge. The fact that strong deviations 
are found only at the highest temperature can be re- 
lated to the existence of two mechanisms controlling j c 
with different temperature dependences. If so, only the 
mechanism dominant at high T has to produce an in- 
homogeneous j c . An example of two such mechanisms 
can be bulk and inter-grain pinning, which are known 
to have different T-dependences£3cJ In thin YBa2Cu3- 
07_5 films the second mechanism can be realized on any 
planar defect such as a boundary between mictoblocks 
with slightly different crystal axis orientationpj a twin 
boundary, or a microcrack.ES 



V. CONCLUSIONS 

The flux density distribution in a superconducting thin 
strip with a general j c (B) is shown to have a special kind 
of symmetry in the remanent state after large applied 
field. Probing the symmetry of measured flux distribu- 
tions is suggested as a simple method to test applicability 
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of the critical-state model without a priori knowledge of 
jc(B). The procedure is simpler than calculation of the 
current distributions because it requires knowledge only 
of the field inside the strip and it is also weakly sensitive 
to "noise" in the experimental data. The method has 
been applied to a thin YBa2Cu307_^ strip which exhib- 
ited a fairly good CSM behavior well below T c , but large 
deviations from the symmetry were observed at 82 K. 
Our analysis shows that the deviations can be attributed 
to an explicit coordinate dependence of j c since both a 
surface barrier, and strong flux creep would break the 
symmetry in a different way. 
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